Novel quinazolines bearing 1,3,4-thiadiazole-aryl urea derivative as anticancer agents: design, synthesis, molecular docking, DFT and bioactivity evaluations

A novel series of 1-(5-((6-nitroquinazoline-4-yl)thio)-1,3,4-thiadiazol-2-yl)-3-phenylurea derivatives 8 were designed and synthesized to evaluate their cytotoxic potencies. The structures of these obtained compounds were thoroughly characterized by IR, 1H, and 13C NMR, MASS spectroscopy and elemental analysis methods. Additionally, their in vitro anticancer activities were investigated using the MTT assay against A549 (human lung cancer), MDA-MB231 (human triple-negative breast cancer), and MCF7 (human hormone-dependent breast cancer). Etoposide was used as a reference marketed drug for comparison. Among the compounds tested, compounds 8b and 8c demonstrated acceptable antiproliferative activity, particularly against MCF7 cells. Considering the potential VEGFR-2 inhibitor potency of these compounds, a molecular docking study was performed for the most potent compound, 8c, to determine its probable interactions. Furthermore, computational investigations, including molecular dynamics, frontier molecular orbital analysis, Fukui reactivity descriptor, electrostatic potential surface, and in silico ADME evaluation for all compounds were performed to illustrate the structure–activity relationship (SAR). Supplementary Information The online version contains supplementary material available at 10.1186/s13065-024-01119-0.


Introduction
Cancer is a complex and severe human disease characterized by abnormal cell growth that can potentially spread to other parts of the body.It is currently recognized as the second leading cause of death worldwide [1].Developed countries are experiencing an alarming increase in cancer-related mortality rates, with approximately 10 million deaths recorded in 2020.Predictions suggest that this number will rise to 16.3 million by 2040 [2,3].
Furthermore, several commercially available anticancer drugs, including Erlotinib, Lapatinib, Vandetanib, Afatinib, and Gefitinib, contain quinazoline as a key component and are commonly used worldwide.
Thiadiazole is another common 5-membered heterocyclic system containing one sulfur atom and two nitrogen atoms.This heterocyclic core structure is present in numerous commercially available drugs and biologically active natural products.The remarkable aromaticity, non-toxicity, and exceptional in vivo stability of thiadiazole make it highly desirable for use in bioactive compounds [22].
The urea functionality plays a crucial role in drug development and medicinal chemistry due to its ability to form multiple stable hydrogen bonds with various protein and receptor targets.Urea and its derivatives have been utilized to modulate drug potency, selectivity, and improve drug properties in the development of anticancer, antibacterial, anticonvulsive, anti-HIV, and antidiabetic agents [32][33][34].
Notably, several FDA-approved drugs containing the urea functionality, such as Regorafenib, Sorafenib, and Tivozanib, have been successfully developed and are widely used as effective anticancer drugs in the global market [35][36][37] (Fig. 1) Considerable attention has been focused on ongoing research dedicated to discovering innovative and effective anti-cancer compounds that can target multiple pathways involved in the growth and spread of cancer cells.The process of angiogenesis, which encompasses the formation and maintenance of blood vessels, plays a vital role in cell development and replication.In cancer cells, this physiological process takes on a similar role and becomes crucial in the progression and metastasis of tumors.Therefore, utilizing efficient anti-angiogenesis agents may be proposed as a suitable approach to combat cancer [38].
Vascular endothelial growth factors (VEGF) and their receptors (VEGFR) have been proven to play significant regulatory roles in pathological angiogenesis, particularly in cancer.Among the three subtypes of VEGFRs, VEGFR-2 plays a critical role in promoting tumor angiogenesis.The interaction between VEGF and VEGFR-2 results in dimerization and subsequent auto-phosphorylation of tyrosine residues (Tyr1059 and Tyr1054) within the receptor.This activation initiates signaling pathways that contribute to angiogenesis in the tumor microenvironment [39].
Hence, the use of antitumor substances that disrupt tyrosine phosphorylation demonstrates inhibitory effects on angiogenesis [40].Therefore, numerous compounds bearing valuable scaffolds particularly quinazolines, thiadiazole, and urea functionality showing great VEGFR-2 inhibitory potencies have been developed.For example, the compound SKLB1002 is a type of quinazoline-bearing 1,3,4-thiadiazole that has been identified as a potent inhibitor of VEGFR-2.
It has shown minimal toxicity in studies [41,42].Diaryl ureas have traditionally been used to inhibit VEGFR-2 activity, but there have been advancements in modifying their structure to enhance their effectiveness.For instance, researchers have substituted the urea functional group with a thiourea moiety and incorporated additional rings such as 1,3,4-thiadiazole, oxadiazole, and 1,2,3-triazole to improve binding affinity with the receptor [22,43,44].
In the present study, the researchers aimed to develop novel quinazoline conjugates with 1,3,4-thiadiazole and diaryl urea, and to investigate their anticancer activities [45][46][47][48][49][50][51].Therefore, an efficient, multi-step, synthetic pathway was executed to produce the desired compounds, which were then assessed for their ability to inhibit cell proliferation in A549 cells (human lung cancer), MDA-MB-231 (human triple-negative breast cancer), and MCF7 (human hormone-dependent breast cancer) in comparison with Etoposide as the reference drug.Moreover, comprehensive computational investigations, including molecular dynamics, molecular docking studies, frontier molecular orbital analysis using DFT calculations, and in silico ADME studies, were performed to confirm the biological results (Fig. 2) To reveal the significant role of the nitro group on this scaffold, it must be noted this functionality is a unique functional group in medicinal chemistry.It has a powerful electron-attracting potency that creates localized electron-deficient locations within molecules and interacts with biological nucleophiles present in living systems, such as amino acids, proteins, nucleic acids, and enzymes.The interaction can happen by nucleophilic addition and electron transfer involving reduction and oxidation, or also simply by molecular complexation, to Fig. 1 Structures of FDA-approved EGFR and VEGFR-2 inhibitor drugs featuring quinazoline, thiadiazole, and urea moieties induce undesired or desired biological changes [52,53].Therefore, countless medicinal chemistry campaigns have been initiated to investigate compounds containing nitro groups.Already, drugs bearing nitro groups have a long history of use as antibiotic, antineoplastic, and antiparasitic agents, as well as fungicides, tranquilizers, insecticides, and herbicides [53][54][55][56][57][58][59].Moreover, some medicinal compounds including PR-104, Tarloxotinib (TH-4000), and Paclitaxel prodrug possess nitro group [60,61].

Materials
The solvents and reagents utilized in this study were obtained from reputable companies such as Merck or Sigma, and no additional purification was required.Melting points were measured using a Kofler hot-plate microscope apparatus, and the values reported were unadjusted.The NMR spectra were obtained using a Bruker FT-500 and Bruker FT-300 MHz spectrometer, with either DMSO-d6 or CDCl 3 serving as the solvent, while TMS was used as the internal standard.The chemical shifts and coupling constants are presented in δ (ppm) and J (Hz), respectively.Additionally, all reactions underwent monitoring via TLC on plastic sheets coated with silica gel 60 F254.The Perkin-Elmer Spectrum Version 10.03.06 was used to obtain IR spectra, specifically using potassium bromide disks.Mass spectra were obtained using the HP Agilent Technologies 5937 instrument with an ionization potential of 70 eV.Elemental analysis for carbon (C), hydrogen (H), and nitrogen (N) was conducted using the elemental analyzer GmbH VarioEL.

Synthesis of quinazoline-4(3H)-one (2)
A mixture of 2-aminobenzoic acid 1 (0.137 g, 1 mmol) and formamide (0.045 g, 14 mmol) was heated with vigorous stirring at 150 ℃ for 6h.After completion of the reaction according to the TLC analysis, the reaction mixture was cooled down to the ambient temperature.The precipitated product was filtered and washed with sufficient amount of water to remove the excess of formamide and to afford pure compound as a white solid with m.p. 212-214 ℃ [62].

Synthesis of 6-nitroquinazolin-4(3H)-one (3)
Under the ice bath conditions, quinazoline-4(3H)-one (0.146 g, 1 mmol) 2 was slowly added into an acid mixture (concentrated sulfuric acid and concentrated nitric acid with a ratio of 1:1) (2 mL) within almost 1h.When the addition was complete, the temperature of the system was slowly raised up to 95 °C, and at this temperature, the reaction was heated for 1 h.After completion of the reaction according to the TLC analysis, the reaction mixture was poured into ice water (25 mL), and it was stirred till yellow solid was precipitated.The precipitation was filtered to afford the pure 6-nitroquinazolin-4(3H)-one 3 with m.p. 279-283 ℃ [63]. 1

Synthesis of 5-amino-1,3,4-thiadiazole-2-thiol (6)
In a reaction flask, thiosemicarbazide 5 (1 mmol) and Na 2 CO 3 (1 mmol) were combined with absolute ethanol (5 mL) and heated to 60 °C for 30 min.A solution of carbon disulfide (CS 2 , 3 mmol) in absolute EtOH (5 mL) was then slowly added drop by drop to the mixture.The resulting mixture was refluxed overnight.Once the reaction was complete, as confirmed by TLC analysis, the solvent was evaporated under reduced pressure.The remaining residue was diluted with water (25 mL), followed by the slow addition of concentrated HCl solution (5 mL).This resulted in the precipitation of compound 6 as a light yellow solid with a yield of 78% and a melting point range of 233-235 [65].

In vitro cytotoxicity assay
The MTT assay was used to test the cytotoxicity of novel 6-nitroquinazoline conjugated with 1,3,4-thiadiazole and diaryl-urea 8 against the A-549 (human lung cancer), MCF-7 (human hormone dependent breast cancer), MDA-MB-231 (human triple negative breast cancer), and HDF (Human Dermal Fibroblasts) cell lines.All of the cell lines were obtained from the Pastor Institute's National Cell Bank in Tehran, Iran, and were cultured in either RPMI-1640 or DMEM supplemented with 10% fetal bovine serum (Gibco, Milano, Italy).The cancerous cells were cultured at 37 °C in a humidified 5% CO 2 incubator after being seeded onto 96-well micro-plates at a density of 1 × 10 4 cells/mL.After overnight incubation at 37 °C, compounds were added to the cells at varying final concentrations (1 µg/mL, 5 µg/mL, 20 µg/mL, and 40 µg/ mL) after being dissolved in dimethyl sulfoxide (DMSO).Fresh phenol red free RPMI-1640 containing 0.5 mg/ mL of MTT was added after incubation.Using a multiwell plate reader (Gen5, Epoch, BioTek) optical density was measured at 492 nm.Both Etoposide and 1% DMSO were employed as positive and negative controls, respectively.Three separate dose-response curves were used to get the IC 50 values (the dosage needed to suppress cell growth by 50%).

Chemistry
The synthetic approach toward the desired novel 6-nitroquinazoline conjugated with 1,3,4-thiadiazole and diaryl-urea 8a-l is illustrated in Scheme 1. Heating a mixture of 2-aminobenzoic acid and formamidine under the reflux conditions led to obtain quinazoline-4(3H)-one 1 which subsequently underwent the electrophilic substitution with the mixture of nitric acid and sulfuric acid at the ambient temperature to give 6-nitroquinazolin-4(3H)-one 2 The chlorination of carbonyl group in this compound with thionyl chloride in DMF was carried out the reflux conditions to obtain 4-chloro-6-nitroquinazolin-4(3H)-one 3. On the other hand, the production of 1-(5-mercapto-1,3,4-thiadiazol-2-yl)-3-arylurea 7a-l was achieved through a condensation reaction involving thiosemicarbazide 4 and carbon disulfide.This reaction took place in the presence of sodium carbonate under reflux conditions in absolute EtOH (compound 5).Subsequently, various substituted phenyl isocyanates 6a-l were added to the mixture in CH 3 CN to yield the corresponding derivatives 7a-l.Finally, chlorine moiety compound 3 underwent nucleophilic aromatic substitution with synthesized 1-(5-mercapto-1,3,4-thiadiazol-2-yl)-3-arylurea 7a-l in the presence of Et 3 N in CH 3 CN under the reflux conditions to afford our desirable 6-nitroquinazoline conjugated with 1,3,4-thiadiazole and diaryl-urea derivatives 8a-l.The spectra of the compounds are provided in the supplementary information file.
The generality of the substrate scope was investigated using various substituted phenyl isocyanates 6a-l bearing electron-donating groups like methyl and methoxy, as well as electron-withdrawing groups like fluorine, chlorine, bromine, and trifluoromethyl.The structures of isolated 6-nitro quinazoline conjugated with 1,3,4-thiadiazole and diaryl-urea derivatives 8a-l are summarized in Table 1.They are deduced based on their IR, 1 H, and 13 C NMR, as well as MASS spectroscopy and elemental analysis.Partial assignments of these resonances are given in the Experimental Part.

In Vitro antiproliferative activity
The targeted 6-nitroquinazoline conjugated with 1,3,4-thiadiazole and diaryl-urea 8a-l were evaluated for their in vitro antitumor activities against three human cancerous cell lines, including A549, MDA-MB231, and MCF7 cell lines by using the MTT colorimetric assay to investigate the role of substituents on the phenyl moiety originated from arylisocyanates.The cytotoxicity is expressed as the concentration which inhibits 50% of cell viability (IC 50 ).In present study, Etoposide was used as a positive control having IC 50 values of 25.8 ± 0.0008 µM against A549, 27.9 ± 0.002 µM against MDA-MB231, and 29.8 ± 0.0007 µM against MCF7.These compounds were shown better cytotoxicity comparable to Etoposide on MCF7 cells, as some derivatives named 8b, 8c, and 8g exhibited IC 50 values of 29.0 ± 0.005 µM, 25.0 ± 0.001 µM, 39.1 ± 0.001 µM, respectively.Other results were noticeably weaker than standard drug.It seems that the presence of methyl as a mild electron-donating group particularly at C-4 position (compound 8c) improved the cytotoxic potency (Table 2).Furthermore, we assessed the cytotoxic impact of 8b and 8c as most potent cytotoxic compounds on normal HDF.Remarkably, no toxic properties were observed in HDF cells even at a concentration as high as 50 µM.

Computational details
To investigate the electronic features of synthesized 6-nitroquinazoline conjugated with 1,3,4-thiadiazole and diaryl-urea 8a-l, quantum calculations were performed at the B3LYP/6-31 g (d, p) level of theory.Frequency calculations were performed to confirm the nature of minima  structures.All the DFT quantum calculations were conducted using GAMESS [67].The frontier molecular orbital analysis was performed to study the chemical reactivity of ligands 8.The electrostatic surface potential (ESP) analysis was performed by Multiwfn 3.3.8code [68].To reveal the relationship between the structure and activity of compounds, the quantitative ESP analysis was done.To evaluate the plausible interactions between most potent compounds, 8c, as a potential VEGFR-2 inhibitor, and VEGFR-2 target, the molecular docking calculation were performed using Autodock 4.2.1.The crystallographic structure of 3WZE was selected as a VEGFR-2 target.Molecular dynamic simulation was performed with GROMACS 2020.Amber force field was used for building the topological files for protein and ligand.To consider the solvation, the TIP3 water model was employed.The box system including protein, ligand, water, and ions was minimized foe 2175 steps.The equilibration process was done in the NVT and NPT ensembles at the pressure of 1 atm and temperature of 300K for 200 ps.The MD simulation was performed for 100 ns.The RMSD, RMSF, and hydrogen bond analysis were performed on ythe trajectory file.

Frontier molecular orbital analysis
The frontier molecular orbitals could be used to elucidate the chemical reactivity of a molecule.The energy gap can be considered as an indicator of the stability and chemical reactivity.Molecules with smaller band gap categorized as more polarizable molecules with higher chemical reactivity.In all 6-nitroquinazoline conjugated with 1,3,4-thiadiazole and diaryl-urea 8a-l, HOMO is located on the phenyl urea moieties whereas the LUMO is localized on the nitro quinazoline moiety.Therefore, according to different distribution of HOMO and LUMO, one can expect that there is a favorable intramolecular charge transfer for these series of ligands (Fig. 3).
The HOMO and LUMO orbital energies can be used to assess the quantum chemical reactivity indexes including electrophilicity and nucleophilicity indexes, chemical potential, softness and hardness parameters as: The quantum chemical indexes including chemical potential, chemical hardness, and global electrophilicity indexes are calculated for all obtained compounds 8a-l, and the results are summarized in Table 3.These parameters are calculated based on the ground state optimized geometries.It can be seen that the chemical reactivity of all ligands is negative which can be assign to the stability of these compounds 8.

Fukuie reactivity descriptor
A key parameter to describe the reactive area in a molecule is Fukui function, f(r), which is based on the conceptual density functional theory.It can be described as where N and υ can be defined as number of electrons in the system and the external potential due to nuclear charges, respectively.Finding the reactive centers is important to elucidate the interaction between the compound and biological target.For this purpose, the Fukui parameters are calculated for the most potent 6-nitroquinazoline conjugated with 1,3,4-thiadiazole and diarylurea, 8c.In general, the centers with high activity show the more values of Fukui functions.This function can be evaluated for three conditions including nucleophilic attack f + (r) , electrophilic attack f − (r) , and radical attack f 0 (r) according to: Dual descriptor is another parameter for the evaluation of the nucleophilic and electrophilic centers at the same time.The regions with �f > 0 and �f < 0 can be described as suitable centers for a nucleophilic attack and electrophilic attack, respectively.f can be evaluated as The graphical view of these indexes was indicated in Fig. 4.
Table 4 The values of Fukui functions f − (r), f + (r), and f 0 (r)) and dual descriptors (∆f(r)) for active sites on the potent 8c ligand Atom q(N) q(N + 1) q(N- Furthermore, to quantify the amount of reactivity in each site of the molecules, the condensed Fukui functions are calculated and the results are reported in Table 4. It can be seen that regarding the nucleophilic attack, the potential active sites are distributed mainly on the nitro quinazoline part.As to electrophilic attack, the  plausible active sites are localized on the phenyl urea moiety.It can be seen that the thiadiazole ring is common in both nucleophilic and electrophilic attacks.To clarify and asses the nature of this position, the dual descriptor parameter could be useful.It can be seen that the DD parameter at thiadiazole ring is mainly positive.Therefore, it has potential for nucleophilic attack.For the other fragments in the compound 8c, the results of dual parameter are in the line with the results of f + (r) and f − (r) Fukui functions.For radical attack, the plausible active sites are located mainly on the nitro and phenyl ring.

Electrostatic potential surface
The electrostatic potential iso-surfaces can be considered as a valuable tool to conduct a close link between the distribution of electrons in a molecule and the biological activity.It can be used to predict and interpret the weak non-covalent interactions between the small ligands and the macromolecule target.The ESP analysis elucidate for us the regions with positive and negative electrostatic potential with low and high amount of electron densities, respectively.The map of electron density for two compounds is plotted in Fig. 5.The regions with low and high electron densities are shown with red and blue color, respectively.The quantitative electrostatic potential analysis was performed to shed light more accurately on the regions with high and low electron densities based on the optimized structures of molecules on the ground state.The positive regions on the hydrogen atoms are due to the lower electronegativities of this element.The electronegative elements with high electronegativities show negative regions.The white color shows region with moderate electronegativities.It is expected that the regions with negative potential conduct strong attraction with positive regions on the target and the positive regions on the ligand conduct attraction interactions with negative regions on the target.It can be seen that different position of methyl group on compounds 8b and 8c led to different distribution of electron density that has an important role in the bioactivity of these compounds.
To have a favorable interaction between ligand and bio target, a complementary electron density distribution must be existed.Different substituents on the phenyl ring led to different distribution of electron density and consequently the different bioactivity.In all 6-nitroquinazoline conjugated with 1,3,4-thiadiazole and diaryl-urea 8a-l, hydrogen atoms with lower electronegativity show positive potential region (red color).The atoms with more electronegativity including O, F, and N show negative potential regions with blue color.

In silico ADME evaluation
The physiochemical and pharmacokinetic properties for all 6-nitroquinazoline conjugated with 1,3,4-thiadiazole and diaryl-urea 8a-l were computed using Swiss ADME online (http:// www.swiss adme.ch/ index.php) toolkit [69].Through this in silico study under the Lipniski's rule, several determined drug-likeness parameters were compared with the known drugs [70][71][72][73][74].These evaluated parameters are summarized in Table 5.In terms of drug likeness properties, it can be seen that there is not any violation to Lipinski rules for these compounds.Consensus Log P o/w is predicted in the range of 2.23-3.85.The calculations predict no Blood Brain Barrier (BBB).

Molecular docking studies
To investigate the interactions between 8c and VEGFR2 target, the molecular docking calculation was performed.There are two types of interactions between SM1 ligand and target VEGFR2.The ligand forms hydrophobic interactions with LEU 840, LEU 889, ILE 892, VAL 899, LEU 1019, LEU 1035, PHE 1047 and hydrogen bonds with LYS 868, GLU 885, CYS 919, and ASP 1046 (Fig. 6).
To assess the stability of the complex including protein and ligand, an MD simulation was performed for 100 ns.It can be seen that, in Fig. 7, the number of hydrogen bonds between protein and ligand vary in the range of 1-4 in most of the time of the simulation.The RMSD plot which is an indicator of stability of the protein, remains stable 0.15 nm after 50 ns.RMSF, which is used to evaluate the mobility of the residues, show low fluctuations and the structure remains stable during simulation time, the peaks are related to the residues out of the active site of protein.Thus, according to the results of RMSD, RMSF, and H-bond interactions, we can see that the structure remains stable during the 100 ns.

Conclusion
A novel series of 1-(5-((6-nitroquinazoline-4-yl) thio)-1,3,4-thiadiazol-2-yl)-3-phenylurea derivatives 8al were designed, synthesized, investigated for their potential as an anticancer agent against A549, MDA-MB231, and MCF7.Among them, compound 8b and 8c bearing methyl group at C-3 and C-4 positions of phenyl ring showed good in vitro anti-proliferative potency with IC 50 value of 29.0 ± 0.005 µM and 25.0 ± 0.001 µM against MCF7, respectively, which was comparable with Etoposide (IC 50 = 29.8 ± 0.0007 µM).The ADME properties were calculated to assess the drug-like features.The quantitative electrostatic potential analysis was done to have a better understanding of the relationship between the electronical structure and activity.The stability of the protein-ligand complex was verified through 100 ns MD simulation.Finally, considering better cytotoxicity of compound 8c, the molecular docking study for this ligand was performed, showing several hydrogen bonds, hydrophobic interactions, and van der Waals forces between it and residues in the active sites of VEGFR-2.Although it has been anticipated that these compounds possess great anti-proliferative potencies, this goal was not achieved which might be related to the adverse role of the nitro group.Therefore, changing this functionality in an attempt to find further potent anticancer agents could be an interesting goal for future studies.

Fig. 4
Fig. 4 Fukui functions for a nucleophilic attack f + (r) b electrophilic attack f − (r) c dual descriptor (Δf), and d radical attack f 0 for 8c ligand at B3LYP/6-31 g (d, p) level.The green and blue regions show positive and negative values of these functions

Fig. 5
Fig. 5 The quantitative ESP analysis on the molecular vdW surface of studied ligands.The unit is in kcal/mol.The green and magenta spheres represent surface local minima and maxima of ESP values, respectively

Fig. 6 .
Fig. 6. 3dimensional (3D) binding interactions of compound 8c.The hydrophobic interactions and hydrogen bonds are indicated as black dotted and blue lines, respectively

Table 5
Drug likeness properties of ligands